On Ferromagnetism in the Large-[7 Hubbard Model 
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We study the Hubbard model on a hypercubic lattice with regard to the possibility of itinerant 
ferromagnetism. The Dynamical Mean Field theory is used to map the lattice model on an effective 
local problem, which is treated with help of the Non Crossing Approximation. By investigating 
spin dependent one-particle Green's functions and the magnetic susceptibility, a region with nonva- 
nishing ferromagnetic polarization is found in the limit U — + oo. The 5-T-phase diagram as well as 
thermodynamic quantities are discussed. The dependence of the Curie temperature on the Coulomb 
interaction and the competition between ferromagnetism and antiferromagnetism are studied in the 
large U limit of the Hubbard model. 



PACS: 75.10.Lp, 71.27+a 

The microscopic description of ferromagnetism in 
narrow-band metals like Fe, Ni, Co and others is one 
of the most interesting problems in solid state physics. 
Since the electrons in these systems are mobile one can- 
not use localized-spin models with effective interactions 
like e.g. the Heisenberg model, but has to take into ac- 
count this itineracy together with the electron-electron 
interaction on a more fundamental level. The first model 
set up to describe such a system is the Hubbard model 



H = -t 



E 

<ij>T 



C icr C jc 



U 



n iT n u . 



(1) 



However, it was realized relatively early that the Hub- 
bard model rather seems to be a generic model for antifer- 
romagnetism and a correlation driven metal-insulator 
transition instead. Just these properties made it an early 
and rather successful candidate for the description of the 
high-T c compounds ||. 

Nevertheless the question about ferromagnetism in the 
Hubbard model was never abandoned, since one of the 
few rigorous theorems about this model definitely proves 
its existence. In 1965 Nagaoka showed that for U = oo 
and one hole doped into the half-filled band the state 
with a fully polarized background is the ground state for 
several lattice structures due to a gain of kinetic energy 
for the hole Q . This theorem initiated a large amount of 
work on questions like the stability of the Nagaoka state 
with respect to doping 8, finite U, etc. |5). Moreover, 
even after 30 years of research the situation appears to 
be rather controversial, especially for bipartite lattices: 
One obtains critical dopings in the range S c from to 0.3, 
depending on the method used . Only for the infinite 
dimensional hypercubic lattice the situation seems to be 
clear: The work by Fazekas et al. suggests that the Na- 
gaoka state is unstable for any finite doping unless 
explicitely favoured by long-range Coulomb interactions 
|jj"of or band structure effects [^l) . 

Most of the above studies are based on a variational 
ansatz and are thus restricted to the Nagaoka - i.e. fully 



polarized - state and its stability as the ground state. 
There is still the possibility of partially polarized ferro- 
magnetism and in any case the necessity to calculate T c 
as function of 5 and U etc. Generally speaking the ques- 
tion to what extent ferromagnetism is a generic feature 
of the Hubbard model or not is still unanswered. 

In this letter we discuss the magnetic phase diagram 
of the Hubbard model on a hypercubic lattice for large 
Coulomb repulsion U. In the latter limit the ground-state 
and low-energy properties of the model are well captured 
by a t-J model with an effective antiferromagnetic ex- 
change J = 2t 2 /U flf]. To solve the t-J model or, more 
precisely, the underlying Hubbard model at U — oo we 
use the dynamical mean field theory (DMFT) jL3|. This 
theory leads to purely local dynamical renormalizations 
of one-particle properties, which can be obtained from 
an effective impurity problem coupled to a self-consistent 
medium 

In addition to the one-particle properties the DMFT 
also allows to calculate two-particle correlation functions 
and thermodynamic quantities consistently |16|]17| ]. Es- 
pecially in the limit of large U one obtains ]17[ 



and 



2t 2 

X U F {T) = X?{T)l[\ + 2d-— xf(T) 



( 9t 2 
X U af(T) = Xaf(T)/ f 1 - 2d- —Xaf(T) 



(2) 



(3) 



for the homogenous (Xf(T)) and staggered (XafCO) 
susceptiblities of the Hubbard model, respectively. The 
quantity Xo'(T) denotes the susceptibility for U = oo 
and d is the spatial dimensions of the system. These ex- 
pressions allow to discuss the influence of finite U once 
the x'S'iT) are known. 

In the follwing we use 4dt 2 = 1 as energy unit. With 
this choice the bare density of states for the hypercu- 
bic lattice is of Gaussian form for large d: po(e) = 
l/v^r exp(— e 2 ). The effective impurity problem of the 
DMFT is solved within the NCA 0,111 and for U = oo 



1 



we furthermore do the calculations with spin-dependent 
quantities to explicitely look at the properties in the 
symmetry-broken phase. An extension of these calcu- 
lations to finite U is extremely tedious and studies along 
this line are in progress ]2Q|. 
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FIG. 1. Magnetization m(T) (crosses), inverse homogenous 
susceptibility Xf (0 _1 (circles) and staggered susceptibility 
XafCO (squares) for 8 — 0.03 as function of temperature. 
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FIG. 2. Magnetic phase diagram for the Hubbard model 
on a hypercubic lattice. The dashed line represents a fit to 
the last data points and predicts a critical doping 5 C « 0.33 
beyond which Tc = 0. 



We begin the discussion of our results with the case 
U = oo. Figure [l] shows the inverse susceptibility (ho- 
mogenous and staggered) as function of temperature T 



for a doping 6 = 1 



0.03. While Xaf( t ) re- 



mains finite for all T, Xf CO ~ 1 vanishes linearly for a 
Tc > and below Tc we observe a finite magnetization 
m = n-f — ri| with m(T) = c^J\T — T c \. Note that the 
critical points found from m(T) and Xj?CO conincide 
indicating a second order transition (see Fig. [l]). Un- 
fortunately our data are not sufficient to extrapolate for 
m(T = 0), which according to the results by Fazekas et 



al. H we would expect to have a value m(T = 0) < n. 
Clearly, this point needs further investigation. 

Repeating the above calculation for different dopings 
we obtain the <5-T-phase diagram in Fig. |^, which shows 
a fairly extended region of ferromagnetism with a max- 
imum in T c at 6 between 0.07 and 0.08. Although the 
NCA in principle does not allow to do calculations down 
to T c beyond 6 = 0.1, the observed Curie- Weiss form of 
Xl?CO enables us to obtain data points in this region 
of the phase diagram from xIf CO a ^ high temperatures. 
Obvioulsy this procedure becomes less accurate for in- 
creasing doping so that the behaviour of the phaseline 
T c {5) currently remains unknown for 5 > 0.2. The ex- 
trapolation of the available data nevertheless indicates 
that a critical doping S c between 0.3 and 0.4 exists be- 
yond which Tc = 0. 

The stability of the ordered phase depends on the inter- 
play of internal energy E(T) and entropy S(T) entering 
the free energy F = E — TS. In Fig. || we thus show the 
difference in free energies AF(T) = F FM {T) - F PM (T) 
together with the internal energy, specific heat and en- 
tropy for the ferromagnetic and paramagnetic state. To 
relate the data to the preceeding discussion the magne- 
tization and inverse susceptibility from Fig. [l] are shown 
again in the upper part of Fig. [| Below Tc the difference 
AF(T) becomes negative, i.e. the ferromagnetic state is 
indeed thermodynamically stable. 
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FIG. 3. Squared magnetization m(T) 2 , inverse homoge- 
nous susceptibility x'FiT) -1 , difference in free energies 
AF(T) = F FM (T) - Fpm(T), internal energy E(T), specific 
heat C(T) and entropy S(T) for 8 = 0.03 in the paramagnetic 
(circles) and ferromagnetic phase (crosses) close to Tc- From 
AF(T < Tc) < it is clear that below T c the ferromagnetic 
solution is stable. 
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The internal energy E(T, n) is given by the expectation 
value of H, which for U = oo is equivalent to the kinetic 
energy. As is evident from Fig. || E(T) for the ferromag- 
netic solution is lower than the corresponding values in 
the paramagnet below T c . Therefore the transition to the 
ferromagnetic phase is obviously connected with a gain 
in kinetic energy. This leads to the conjecture that the 
physics underlying the stability of the ferromagnetic state 
should be roughly the same as in the particular case stud- 
ied within the Nagaoka theorem. Since at T c the slope of 
the internal energy E(T) changes for the ferromagnetic 
solution the specific heat C(T) = dE/dT shows a jump 
characteristic for a second order phase transition. Note 
that only a very small temperature region around Tc is 
shown in Fig. y and that C(T) decreases again for lower 
temperatures. Finally, the entropy S(T) is obtained from 
F = E — TS. Just above Tc its value is very close to 
In 2, the value expected for a spin 1/2 system, while below 
Tc the increasing spin order leads to a strong decrease 
of S(T). 




FIG. 4. Spectral functions for 5 = 0.03 and (3 = 70 for 
both spins. The temperature is well below T c , so that the 
spin dependent solution (dashed line) shows a difference in 
the spectral weight for the two spin directions. The inset 
shows details near the Fermi energy. 

The ordered state of course also shows up in the dy- 
namical properties such as the one-particle Green's func- 
tion. In Fig. H we show the density of states (DOS) for 
a doping S = 0.03 at a temperature T = 1/70 < Tc 
for both the paramagnetic (full line) and ferromagnetic 
solution (dashed lines). The basic features in the ferro- 
magnetic phase are similar to those of the paramagnet. 
One finds the lower Hubbard band represented by a broad 
peak and a quasi-particle resonance near the Fermi en- 
Due to U 



ergy 14 lq|. Due to U — oo the upper Hubbard band 
does not appear. 

Depending on temperature and doping spectral weight 
is transferred between the states with spin a and that 
of — a, most prominent in the charge-fluctuation peaks, 



which results in differences in the occupation numbers n a 
and in a finite magnetization m = nj — nj,. Note that this 
does not occur due to any explicit magnetic exchange but 
rather to the fact that the energy loss by increasing the 
population of the er-band is outweighed by the gain in 
kinetic energy from the holes in the — cr-states 0. In ad- 
dition the peak positions for the minority/majority spins 
are shifted to a somewhat higher/lower energy In terms 
of a band picture this means a slight splitting of the lower 
Hubbard band. 

The energy splitting is observed for the quasi-particle 
resonance near the Fermi energy as well (see inset of 
Fig. ||). In contrast to the lower Hubbard band both 
peaks show a loss of spectral weight compared to the 
paramagnetic state. This reflects the suppression of the 
Kondo like effect underlying the quasi-particle resonance 
by ferromagnetism, analogous to the effect of an external 
magnetic field in conventional Kondo physics. With in- 
creasing magnetization the resonances will continuoulsy 
decrease in height and eventually vanish for very low tem- 
peratures. 

Let us now turn to the interesting question of the de- 
pendence of these results on U. Generally speaking we 
expect that with decreasing U the Curie temperature Tc 
should be supressed and, due to the antiferromagnetism 
favoured by a finite U, a competition between ferromag- 
netic and antiferromagnetic order should occur. This 
anticipated behaviour can readily be read off the signs 
in equations (||) and (0), i.e. U < oo tends to supress 
X^(T) and enhance x^ F (T). In addition, the Curie- 
Weiss form of Xf{T) and Xaf(T) (see Fig. Q) allows to 
rewrite eqs. (|^) and (^ in such a way as to identify a 
Curie temperature T C (U,S) = Tg>(5) - C{8)/U and a 
Neel temperature T N {U,5) = -Q{5) + C(S)/U, where 
T C °(S) and C(S) are the Curie temperature and Curie 
constant for U — oo, while Q(S) and C(5) denote the 
intercept and inverse slope of XafC^) -1 - A detailed dis- 
cussions of C(6), C(8) and Q(S) will be given elsewhere 
pof . Here we want to focus on the resulting S-U phase 
diagram in Fig. |^, curves A and B, which were obtained 
by plotting max(T c (C/, 5),T N (U, 5),0). 

One sees that for large U an extended region of ferro- 
magnetism exists above curve A, which is completely su- 
pressed for U < U c w 20. For decreasing U the ferromag- 
netic order is eventually replaced by antiferromagnetism 
in the region below curve B. Note that up to 6 as 0.07 
we find a direct transition from the ferromagnet to the 
antiferromagnet, which we would expect to be of first 
order ending in a second order critical point. A more de- 
tailed investigation of the region would thus be of great 
interest. However, since the transition temperatures are 
already very small there we do not see any chance to 
achieve this with the methods currently available. Be- 
yond S > 0.07 a paramagnetic region separates the two 
phases. In addition to our new findings we also include 
results on the phase line between antiferromagnet and 
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paramagnet for the full Hubbard model (1) at small U 
(curve C) |16|. The behaviour for the largest U values 
in this case extrapolates nicely to our phase line B for 
U — > oo. For decreasing U and increasing 8, however, 
the approximation of the Coulomb term by an effective 
exchange becomes worse, i.e. the magnetic order is much 
stronger supressed by doping for a given U. 
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FIG. 5. 5—U phase diagram of the Hubbard model from the 
large U limit (curves A and B). We find a region with ferro- 
magnetism above curve A. Below curve B antiferromagnetism 
is predicted. The phase line at small U (curve C) is taken 
from flu] and extrapolates to our phase line for U — » oo. 



To conclude, we have shown that the Hubbard model 
on a hypercubic lattice provides a scenario of ferromag- 
nctism for finite doping and at finite temperatures. We 
were for the first time able to obtain sensible results for 
Tq as function of 8 and U for the strong coupling case. 
As in Nagaoka's case the phase transition originates from 
a gain of kinetic energy, as we could see from thermo- 
dynamic quantities. The S-T-U phase diagram shows 
a fairly extended region of ferromagnetism for large U 
that is completely suppressed for U < U c ~ 20 and 
8 > 8 C w 0.3. For small doping and large U we observed 
in addition a direct transition from the ferromagnet into 
an antiferromagnet as function of U . 

Unfortunately our method to solve the DMFT does 
not allow to study temperatures T « T c . Thus several 
important questions have to remain unanswered: What 
is the ground-state magnetization m{T — > 0) (cf. j|f§) 
and of what nature is the ferromagnetic «-> antiferromag- 
netic transition, for example. In future work one also 
must investigate the order of the transition paramagnet 
<-> ferromagnet, which is under current discussion (cf. 
ref. |22|| , where a first order transition is stated within a 
different method), more closely. 

We thank Ch. Helm and W. Heindl for helpful dis- 
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